Theory of Superconductivity in the Cuprates 
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The quantum critical fluctuations of the time-reversal breaking order parameter which is observed 
in the pseudogap regime of the Cuprates are shown to couple to the lattice equivalent of the local 
angular momentum of the fermions. Such a coupling favors scattering of fermions through angles 
close to ±7r/2 which is unambiguously shown to promote d-wave pairing. The right order of mag- 
nitude of both Tc and the normalized zero temperature gap A/Tc are calculated using the same 
fluctuations which give the temperature, frequency and momentum dependence of the the anoma- 
lous normal state properties for dopings near the quantum-critical value and with two parameters 
extracted from fit to such experiments. 

PACS numbers: 



I. INTRODUCTION 

The objective of a microscopic theory of the phenom- 
ena in Cuprates ought to be to derive their universal 
properties, in all the parts of their phase diagram, based 
on a single set of ideas and with a consistent set of cal- 
culations on a well-defined model. In particular, since 
superconductivity is an instability of the normal state 
of a metal, it is necessary that the same theory which 
seeks to explain high temperature superconductivity in 
the Cuprates also explain their remarkable normal states. 
To be convincing, the theory should also lead to unique 
predictions which can be tested in experiments. 

Towards these goals, one particular approach to the 
theory has so far achieved the following : 

(1) Starting from the three-orbital model with on-site 
and nearest neighbor repulsions, it was predicted that 
the pseudogap state occurs through a phase transition to 
a new state of matter with spontaneous orbital currents 
without changing translational symmetry. The transition 
temperature T*{x) — > for x — > Xc, defining a quan- 
tum critical point at x = Xc- The loop current order 
parametei— has by now been observed in the pseudogap 
region of four distinct families of CupratesSr^^. 

(2) The quantum critical fluctuations (QCF) of the ob- 
served order have been derived^ and shown to have the 
spectrum with uj/T scaling and spatial locality, which 
was introduced phenomenologically long ago^ to explain 
the observed normal state anomalies^ and predict the 
single-particle spectra in the 'Marginal Fermi-liquid' re- 
gion around x = Xc- Using ARPES data the parameters 
of the QCF's and their coupling strength to fermions have 
also been determinecU^. 

In this paper, we show that the derived QCF couple 
to (the lattice equivalent of) the local orbital angular 
momentum of the fermions. This is a natural general- 
ization of the coupling of the spin-angular momentum of 
the fermions to collective spin-fluctuationsii. We derive 
the momentum-dependence of the coupling of the QCF 
to fermions of the conduction band and show from mi- 
croscopic theory as well as symmetry considerations that 
their exchange leads to an attractive d-wave pairing. Us- 



ing the parameters extracted from the quantitative fit to 
the normal state anomalies and Angle- Resolved Photoe- 
mission Spectra (ARPES), we obtain the right order of 
magnitude of Tc and A for superconductivity in the re- 
gion dominated by the QCF. The principal findings of 
this paper can be tested in detailed analysis of ARPES 
data in the superconducting state using the generaliza- 
tion of the McMillan- Rowell procedure for s-wave super- 
conductors. 



A. Plan of this Paper 

In order to present the new results of this paper, it 
is useful to briefiy recapitulate earlier work upon which 
it builds. In Sec. II, we summarize (i) the microscopic 
model, (ii) the derivation of the Loop-current order based 
on it, and (iii) the quantum statistical mechanical model 
for the quantum fluctuations of the order parameter. As 
shown earlier, lattice anisotropy is irrelevant in the fluc- 
tuation regime and the spectra of the fluctuations is ob- 
tained from the solution of the dissipative quantum xy 
or rotor model in the continuum. In Sec. Ill, we present 
the coupling of the fermions to the fluctuations based on 
general symmetry considerations and show that in the 
continuum, the fluctuations of the angular momentum 
of the rotors couple to the local angular momentum op- 
erator of the fermions. In Sec. IV, we derive, through 
microscopic calculations, the coupling of the fermions to 
the fluctuations in the lattice model and show that it is 
a generalization of the continuum model to take into ac- 
count the lattice anisotropy. Some technical details of 
the lattice calculations are given in the Appendices. In 
Sec. V, we present the vertex for superconductivity de- 
rived from the couplings in Sees. Ill and IV and prove 
that only d-wave pairing is possible in the model. In Sec. 
VI, we derive parameters of the model from the fit to the 
normal state spectral function A(k, lu) of the fermions 
and estimate Tc and the superconducting gap A based 
on generalization of the Eliashberg theory to d-wave su- 
perconductors. We also discuss the limit of validity of 
the Eliashberg theory in the present context. 



II. SUMMARY OF WORK LEADING UP TO 
THE PRESENT WORK 

A. Microscopic Model 



the area formed by the closed loops. Such Gauge in- 
variant combinations of (j}i{R, R') within a unit-cell i are 
organized in to irreducible representations /.t of the point 
group symmetry of the lattice: 
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FIG. 1; The unit-cell chosen for the two-dimensional model 
for Cuprates used for calculations in this paper. The labeling 
used to denote the Cu orbitals and the oxygen orbitals will 
be used throughout the paper. 

The relevant microscopic model for the loop-current 
order (LCO) in the Cuprates is the three orbital modeli^ 
with local and finite range interactions. A unit-cell with 
a Cu and two Oxygen orbitals per unit-cell and with la- 
belling used in this paper is shown in figH] The nearest 
neighbor copper-oxygen interactions, 

Hnn= Y. ypdMR)[npAR') + np^iR') (1) 

{R,R') 

+Vppnp^{R)np^{R')]. 

play a crucial role in deriving the LCO. Here 
nd{R),np^^ {R') are the charge operator for the d-orbital 
on Cu site and the px^y-ovhital on the oxygen sites at 
the sites R and R' respectively. Only neighboring Cu 
and Oxygen sites on the 8 links {R, R') per each unit cell 
are summed. In the metallic state, the local interactions 
(U's) are assumed to only renormalize the kinetic energy 
parameters, unlike in the insulating-AFM state close to 
half-filling. The possible novel changes in symmetry^ are 
seen by re-expressing Eq. ([1]) using the operator identity. 



2aUR)a,{R)bl{R')b^,{R') 

= -|C-(i?,i?')l'+4(-RK(^) 
0^^,{R,R') = iaUR)ba'{R') + h.c. 



bl,{R')b,,iR'), 
(2) 



0{R, R') = J2a ^o-o-(^' R') is proportional to the current 
operator on the link between site i and i' . 

Suppose an expectation value < 0{R, R') >= x{R, R') 
were shown to exist. Given a kinetic energy coefficient 
t{R, R') in the bare Hamiltonian, this is equivalent to an 
effective kinetic energy operator on the link with a com- 
plex coefficient (= t{R,R') + ix{R,R')). This amounts 
to a vector potential 0(-R, R') = arctan(a;/i) on the link 
(-R, i?'). Gauge invariant combinations of vector poten- 
tials always form closed loops and correspond to flux in 



As may be seen from Fig. ([T]) and Fig. Q, there are 5 
closed loops which can be formed in a unit-cell through 
connecting nearest neighbor Cu-O and 0-0 links. This 
is consistent with the vector potentials on each of the 8 
links per unit-cell and the three lattice points per unit- 
cell at which independent gauge transformations can be 
made. The algebraic representation as well as their flux 
patterns of the 5 varieties of the n are exhibited in Ap- 
pendix [Cj In the ground state in the pseudogap phase, 
one of these symmetries (with two-dimensional represen- 
tation Eg) is realized globally in the experiments and also 
found to be the lowest energy state in mean-field or better 
calculationai^. There are differences in the experimental 
results from the prediction of the two-dimensional model. 
However, the spatial symmetry and the number of al- 
lowed configurations per unit-cell in the more general 
three dimensional modeli^ are not changed from those 
given by the two-dimensional modeli^. 

The observed long-range ordered state is described by 
a time-reversal odd polar-vector L which has four possi- 
ble orientations. These four domains are shown in fig ID 
As discussed earlier—, the generalization to the three di- 
mensional solid including apical oxygens does not change 
the symmetry classification or the number of possible do- 
mains of the state. The most important fluctuations in 
the model in the quantum-critical regime are due to the 
transitions between the four configurations possible in 
this regime in each unit-cell. To describe the fluctuations, 
we must in the minimal description retain the term which 
causes transitions among the four local configurations. 





FIG. 2: The four domains of the circulating current phase 
are shown. They may be specified by the four orientations of 
a vector L shown in red. 



Effective Hamiltonian 



III. QUANTUM FLUCTUATIONS 



The complete Hamiltonian obtained using the identi- 
ties in Eq. ([21 O is explicitly derived in appendix (|D|) . 
Limiting ourselves to the most important fluctuations, 
the derived effective Hamiltonian from which earlier re- 
sults and the results in this paper are derived is 



The quantum model is specified in terms of operators 
Li — e*®* , whose eigenstatcs arc the four angles 9i in each 
cell i depicted in Fig. ([2]) : 



jdi 



\0^ 



(8) 



H = K.E. - 
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where 



The labels 



E 



{'-•i,x', '-'i,y' ) 



(4) 



(5) 



will denote the unit-cells. K.E. is the ki- 
netic energy operator which is explicitly written down in 
appendix IBJ in terms of two parameters tpd,tpp, which 
are respectively the effective hopping parameters (renor- 
malized due to on-site repulsions) between the nearest 
neighbor Cu-0 and 0-0 sites respectively. £ stands for 
x',y' = (x±y)/\/2. T-ii.e arid L^^^ are given in terms 
of the fermion operators in the Appendix [C] The oper- 
ators Ijz,i introduce quantum fluctuations in the model 
through transitions between the four possible local con- 
figurations of the local order parameter L, as derived in 
appendix (|E| and discussed below. 

The operators Li, hz^i are expressed as the sum of a col- 
lective part Li,Lz^i and a residual incoherent fcrmionic 
part L/i, hzji using the Hubbard-Stratanovich or equiv- 
alent methods: 



L, 






(6) 



The mean-field Hamiltonian consists of the kinetic en- 
ergy Eq. ljBip and the second term of Eq.Q. Integrat- 
ing over the fermions generated a mean-field free-energy 
F,„/((L))i. Minimizing this gave the stable long range 
order with an order parameter < L^ >=< L^ > for all i. 

In this paper we are concerned with the coupling of 
the quantum critical fluctuations of the order parameter. 
The effective Hamiltonian for such fluctuations is gen- 
erated from eqn. Qj using the substitutions eqnlHl The 
fermion operators in the bilinear terms LiiJjfi + JjjJjfj are 



eliminated by integrating over the propagator ( Lt^L/j 

in the standard manner to generate a coupling between 
the collective variables L^ and L ^ . Keeping coupling only 
between nearest neighbor cells this generates the fluctu- 
ation Hamiltonian, 



{ij) 



(7) 



Here J is in general a second rank tensor whose compo- 
nents depend on the orientations of L^ and Lj. Only 
its order of magnitude can be estimated and is of 
O(10-2y2/*)> where V is Vpd 
nctic energy parameters 
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Vpp and t arc the ki- 



Given the symmetries of the four domains, the angles 
correspond to 



|0 = 7r/4) = 


\x + y) 


5i = 37r/4) = 


\-x + y) 


9 = 5^/4) = 


\-i~y) 


5i = 77r/4) ^ 


\x - y) 



(9) 



The low-energy quantum fluctuations of the order pa- 
rameter form current loops of all sizes and shapes and 
varying in time. They are generated by the elementary 
process of quantum-flips between the four-configurations. 
Wc show explicitly in Appendix lEl that the operator L^.i 
is the generator of rotations in the space of the four one- 
particle eigenstatcs of the operator Li. In other words 
the operator that rotates the states is 



U = exp 



(-'?■ 



1 - iL, 



l!. 



(10) 



The operators U causes transition between \9) and 16* 

^2): 



v,\e,) = \{e-Ti/2),) 



Continuum Model 



(11) 



The functioning of the rotation operator is more fa- 
miliar in the continuum model, where the four states are 
replaced by a continuum of angles 9i. In fact, in the 
fluctuation regime, the discreteness of the Oi variables is 
a (marginally) irrelevant perturbation and a continuous 
distribution of Oi gives the correct correlation functions. 
The model is then just the quantum rotor model : 



H. = Y.^-^ + jY.^os{e.-(^^) 



(12) 



where Lzi = id/d6i, conjugate to the operator 9i, 
causes rotations of 9i. The quantum critical fluctu- 
ations are calculated by supplementing the quantum 
rotor model by the dissipation term of the Caldeira- 
Leggett form. cxp{i9{r,t)) The Fourier transform 
of the spectral function of the correlation function 
(exp(i6'(r, t)) cxp(i6'(r', t'))) derived ini is 



lmxiq,uj) 



-Xotanh(w/2r), |a;| 



0, 



w 



< 
> 



UJr, 



(13) 



The value of the cut-off Wc and of the amphtude xo will 
be deduced from experiments below. As noted this is 
of the same form as suggested phenomenologically. The 
spectral weight xo in that equation may be fixed from 
^ /dw[— Imx(q, w)] « (2$o)^, where $0 is the ordered 
flux in each of the two Cu-0-0 triangular plaquettes in 
each unit-cell. $0 is given from experiments^!^ of an or- 
dered moment of about O.l/is per triangular plaquette. 
In appendix [F1 we show that the singular part of the cor- 
relations of the " angular- momentum" operator, Jj^i are 
proportional to those in Eg. p^ . 



IV. COUPLING OF FLUCTUATIONS TO 
FERMIONS 

A. Coupling in the Continuum Model 

It is instructive to write down the coupling for the 
special case that the conduction electrons are considered 
well approximated by those in the continuum. This can 
be done by symmetry considerations alone; the micro- 
scopic results for the lattice are shown below to reduce 
to these in the continuum limit. 

The minimal coupling of the operator for the angular 
momentum of the collective fluctuations Lj,j to fermions 
can only be to the local angular momentum of contin- 
uum fermions. Thus in the continuum limit, the coupling 
Hamiltonian of the fluctuations to the fermions is 



ex drf{\r\)tJj+{r)irxp)iP{r)L^{r) + h.c. (14) 



Hr. 



where p is the momentum operator, so that (r x p) is 
the angular momentum operator. /(|r|) is a function 
which restricts the integrals to be only over the (circu- 
lar) Wigner-seitz cell of the continuum problem. Fourier 
transforming, we get that 

Hcoup oc 7(k,k')^+(k)V(k')L,(k-k')-fi/.ai5) 
7(k,k') ex i(kx k') 

The important point about Eq. (|15p is that scatter- 
ing of fermions through an angle near 7r/2 or — 7r/2 is 
strongly favored compared to backward tt or forward 
angles. The other important point is the factor of i sig- 
nifying coupling to time-reversal breaking fluctuations. 
These two points are crucial to the the pairing symme- 
try favored, as shown below. It should also be clear that 
what has been derived is effectively the equivalent for 
coupling of collective modes which transform as orbital 
magnetic moments to the orbital moment of fermions, to 
the familiar coupling Jip (r) aip{r) ■ S(r) of coupling be- 
tween collective spin-moment variables S to the fermion 
spins. Note that the physics of d-wave being favored 
through strong enough Antiferromagnetic fluctuations is 
also related to the scattering through an angle near 7r/2 
or — 7r/2 and a minus sign (i^) due to spin-trace of the 
fluctuations of S. 



B. Coupling of fluctuations to Fermions for the 
lattice model 



Now we return to the lattice model of Eq.Q and 
Eq. © to generate the coupling of the collective variables 
to the fermions in microscopic theory. Putting the latter 
into the last term of the former generates the requisite 
coupling term: 



Hr. 



y^ -j^Lz,jjL"'"z,i + h.c. 



(16) 



We can now Fourier transform to get the coupling Hamil- 
tonian in momentum space for states in the conduction 
band. Let us take the simplest representation of the con- 
duction band states (in the absence of orbital order) for 
which the annihilation creation operators are 

Ck.cr = —7={dl^,a + iiSxCk^Pkx.cr + Sy{k)pky,a)s~y) . (17) 

Here Si^(k) = sin(fc2,a/2), Sj,(k) = sm{kya/2), s,j.y(k) = 
Jsl{\i) + Sy(\i.)). Projecting to these states, we 

have 4 ~ 4/^/2, pL ~ ''t{s^{k)/^/2sj;y)cl 
and ply ss -i{sy{'k)/^/2sxy)cl, where s^-y = 
■\/sin(fca;/2)^ + sin(A;j,/2)^. The coupling Hamiltonian is 

Hcoup = J2 7(k,k')4(k')c.(k)*L,,q, (18) 

k,k'CT 

where the coupling matrix is 

7(k,k') = Z (^^) [Sx{k)Syik') - Sy{k)Sx{k')]Sxy{k,k') 

(19) 
Here q = k - k' and S^yik, k') = {s^^{k) + s-^{k')). 

V. PAIRING SYMMETRY: 



Integrating over the fluctuations in Eq. (|T8)) in the stan- 
dard manner gives an effective vertex for scattering of 
fcrmion-pairs: 

Hpa^r^ng- ^ A(k, k')c),, (-k')4 (k')Ca (k)c,. (-k); 
kfrk'cr' 

A(k, k') =7(fc, k'M-k, -fc')Rcx(c^ = fk - 4)- (20) 

The susceptibility appearing in the coupling A(k, k') is 
given in eq[T31 For further discussion see appendix [F] As 
will be discussed, this is correct to 0{Xujc/Ef), where A's 
are dimensionless coupling constants derived below. 

It is illuminating to note first the symmetry of the 
favored pairing due to such a coupling in a continuum 
approximation for fermions near the fermi-energy. In this 
approximation, Sa:(k) (x {kxa)/2, etc. so that 7(k, k') ex 
i(k X k'). The pairing vertex is then 

A (k, k, ) ex -(k X k')^Rex((k - k'), to). (21) 



Since Rex(k— k',a;) < for —lUc < uj < ujc, independent 
of momentum, the pairing symmetry is given simply by 
expressing (k x k')^ in separable form : 

(k X k')' = 1/2 \{kl + kl){k,^ + fc'^) - {kl - kl){k,^ - P) 



4(fc,fc,)(fc;fc;)] . 



(22) 



Pairing interaction in the s-wave channel is repulsive, 
that in the two d-wave channels is equally attractive, 
and in the odd-parity channels is zero. The factor i in 
7(k, k'), present because the coupling is to fluctuations 
of time-reversal odd operators, is crucial in determining 
the sign of the interactions of the pairing vertex. 

We now return to Ea. (j20[) for a semi-quantitative anal- 
ysis of pairing in the lattice. The dimensionless coupling 
constants Aq, determining T^ and the gap A are given in 
a generalized Eliashberg theory byi^ 



A. 



2iV(0) / dMf7'A(k,k')F„(k)i^„(k') 



(23) 



and how ARPES experiments have been used to deter- 
mine the parameters used later in this paper. We argue 
that since superconductivity is an instability of the nor- 
mal state which occurs at T < Tc, it is unlikely that any 
other fluctuations can dominate in determining Tc. 

The value of the cut-off Wc, of the amplitude xo in 
Eq. [13] and the coupling constants in front of 7 in Eq. 
1201 can be provided in terms of the parameters of the 
microscopic model as well as deduced from experiments 
below. Given a coupling function of such fluctuations 
to fermions 7(k, k') to scatter from k to k', calculated 
below, the self-energy of the fermions isi^ 



Iml](w, k) 




(26) 



Here A(k) = A^(0) < 7^ >k' and < 7^ >u' is the average 
of |7(k, k')p over k' on the fermi-surface. In the phe- 
nomenological approach^ this was taken to be momen- 
tum independent. In Ref. ljlOl ). this expression is com- 



Here dO, = N{Q)^^dSk/\vk\ where Sk is an element the 
Fermi surface and A^(0) is an effective density of states on 
the Fermi surface for one spin species. Also, FQ,(k)FQ.(k') 
are separable lattice harmonics on to which A(k, k') is 
projectedi^: 



Fs(k) =A^s[l — cos(fca.a) cos(fcj,a)]; 
i^di(k) =Ndi[cos{kxa) - cos{kya)]; 
-Fd2(k) =A^d2 sm{kxa) sin{kya) . 



(24) 



Where the labels (s, dl, d2) represent the irreducible lat- 
tice representations {Aig,Big,B2g) of a tetragonal lat- 
tice, popularly referred to as extended s wave, d^^-y^^ 
and dxy symmetries, respectively. The factor N^ ensures 
normalization J dilFa{k)Fa{k) = 1. 
The resolution of A(fc, k') in Eal^is 



A(fc,fc') = Ao 



rF,(k)F,(k') F,i(k)F,i(k') F,2)(k)Fd2(k') 



N? 



Nl, 






(25) 

From (US]), the s-wave interactions are repulsive while 
the interaction is equally attractive for the (i(a;^ — V^) 
and dajy-waves for a circular fermi-surface. For the actual 
fermi-surface of the cuprates in which the fermi- velocity 
is largest in the (1, 1) directions and the least in the (1, 0) 
or the Cu-0 bond-directions, d(a;^ — j/^)-pairing is favored 
because in that case the maximum gap is in directions 
where the density of states is largest. 



VI. DEDUCTION OF PARAMETERS FROM 
ARPES EXPERIMENTS AND ESTIMATES OF T^ 

AND A 

In this section, we first summarize experimental evi- 
dence and calculation directly showing that the scatter- 
ing of fermions is uniquely given by the derived QCF's 




Energy (eV) 



FIG. 3: The MDC linewidths along the 7r,pi directions for 
all the measured cuprates. The detailed references for each 
cuprate are given in Ref. jlCI ) 

pared with the data in all available directions and parts 
of the phase diagram of the cuprates. In fig. ([3]), we show 
the deduced MDC linewidth in the (tt, tt) direction for all 
the cuprates near optimal doping for which data is avail- 
able. This data is taken with poor energy resolution, 
> 'iOmeV, to cover a wide energy range. The linearity 
of the linewidth with energy in the normal state for low 
energies has been checked with better precision in other 
experiments. Here we focus on the full energy range. We 
notice the remarkable correspondence with Eq. ()26p with 
the cut-off ujc between 0.4 and 0.5 eV for all the measured 
cuprates. Below we will use the slope of these curves for 
w < Wc to deduce the coupling constant Xi in different 
angular inomentum channels. The normal state resistiv- 
ity and optical conductivity can also be calculated using 
the values of A^ and Wc to within about 30% of those 
deduced from the single-particle spectra. 

The most important point to be noted from fig. ^ 
is the following: The result in Eq. ([25]) arises because 



the scattering at any energy w >> T is proportional 
to the integrated weight of fluctuations up to uj, i.e 
ex /p Inix(a;'). Therefore the linearity of the scattering 
rate with lo up to about lUc and constancy thereafter is a 
direct proof of the fluctuations of the form of Eq. ([T3| . 
The rather sharp ujc proves that one need not be con- 
cerned that a distinct energy scale of fluctuations may 
not existi^. 

An important deduction from recent analysisi^ of high 
resolution laser based ARPES at different angles on the 
fermi-surface is that a momentum independent fluctua- 
tion spectra is obtained from the inversion of the data 
through Eliashberg equations to fit the data at different 
angles. 

In Eq. (P5|) . Aq may be estimated thus: Rex(a;) w 
— 2xo In rj\ for |a;| < uJc and it is vanishing beyond. The 
cut-off Ldc is important but the weak dependence on |w| 
may be ignored by replacing it by 0(Tc) for estimates of 
parameters determining T^. For uJc ~ OAeV, as deduced 
from experimentsiS and Tc ~ lOOK, Rcx(a;) is then sa 
—6x0 for oj < UJc- Using this estimate Aq ~ —6V^N{0)xo 
with a cut-off in the range of interaction at tOc- (Here 
s'^y{k = kp) ~ 1/2 has been used.) 

We now ask whether the Uc and A deduced from exper- 
iments in the normal state yield the right order of mag- 
nitude of both Tc's and the ratio of the zero-temperature 
gap A to Tc- Compared to s-wave superconductors, 
the normal self-energy and inelastic scattering lead to 
a stronger depression in T^, and a stronger enhance- 
ment in A/Tc in d-wave superconductors^^; estimate, 
Tc ~ ajcCxp(— (1 + |As|)/|Ad|), where Ag is the coupHng 
constant which appears in the normal self-energy and A^, 
the coupling constant which appears in the anomalous 
self-energy. From Eq. ^, \>^d\/\Xs\ ~ 1/2. The nor- 
mal state self-energy deduced from experimentsi^ gives 
I As I ~ 1 and uJc ~ OAeV. The formula above gives 
Tc ~ 80K for these values. 

The value of A/Tc may be read off from Fig. (4) 
of Ref. (fish using the parameters above to be about 
2.5Ao/rcOj where Aq and Tco are estimated ignoring the 
self-energy and inelastic scattering. This gives A/Tc ~ 5. 
A/Tc of 4-5 are reported for the cuprates^S. 

We can roughly estimate the value of A's and Wc from 
the microscopic parameters. As estimated above, |As| ~ 
6N{0)V'^Xo- To get A,, w 1 requires Xo « 6 x 10-'^{eV)-^ 
for N{0) sa 2 states /eV/ unit ~ cell,V « leV. For the 
deduced ujc ~ OAeV, this requires an ordered $0 pcr 
triangular plaquettes of about 0(0. 1). As mentioned al- 
ready this corresponds well with the measured moment 
deduced by Fauque et alji^,. 

A comment should be made on the validity of a theory 
of superconductivity with electron-electron interactions 
using an Eliashberg type simplificationiS. The Eliash- 
berg simplification of the theory works only if there is a 
small parameter which limits the irreducible vertex in the 
particle-particle channel to ladder diagrams alone with 
the simplest single-particle self-energy in the propaga- 
tors. For any model with single-particle self-energy which 



is nearly momentum independent in the normal state as 
in the present case (and the case of electron-phonon in- 
teractions), this parameter is 0{Xujc/W), where W is 
the bare electronic band-with. This parameter for the 
cuprates with W « 2eV is 0(l/5). This does not allow 
the luxury of a parameter of O(10~^) as for electron- 
phonon interactions but small enough to have a system- 
atic theory. 

The results of this paper and their applicability for 
a consistent calculation of superconducting parameters, 
can be tested in detail through inversion of Angle Re- 
solved Photoemission (ARPES) data by a procedure, 
which is a generalization of the Rowell-MacMillan pro- 
cedure for s-wave superconductors^!. (Approximate in- 
version of such data has recently appeared^^ with results 
consistent with the results of this paper for the shape of 
the spectrum, and its cut-off energy.) 



VII. CONCLUDING REMARKS 

This paper has derived that the fluctuations responsi- 
ble for the normal state anomalies near the QCP couple 
to fermions to promote d-wave pairing and with parame- 
ters taken from fits to the normal state data and consis- 
tent with estimates from microscopic calculations calcu- 
lated the right magnitude of both Tc and A/Tc. The criti- 
cal fluctuation spectra develops a low-energy cutoff which 
increases as x increases. So Tc is expected to fall. For un- 
derdoping, superconductivity can only be calculated in a 
state with the competing order parameter whose strength 
goes up as x decreases from Xc- Again Tc must fall. The 
details of such calculations arc work for the future. 

Useful discussions with Han- Young Choi and Lijun 
Zhu are gratefully acknowledged. CMV's research was 
partially supported by National Science Foundation 
grant DMR- 0906530 



Appendix A: Introduction to the Appendices 

There are three principal purposes of these appendices. 
First is to derive the Hamiltonian, Eq. ((4]) from the in- 
teractions of Eq.(IT|). Towards this end, we first give the 
expressions for the Kinetic energy so as to specify the 
choice of phases made in the d and the p-orbitals and 
of the current operators in the gauge chosen for the cal- 
culations. This is followed by a detailed derivation of 
Eq.(|3]). The second purpose is to show that the operator 
hiz provides the kinetic energy which acts as a rotation 
operator for the four possible flux patterns given by L^^ 
in the cell i. The third is to show that the correlations 
of the operators L^^ have the same singularities as the 
order parameter fluctuations calculated in Rcfs. (I7ll23h. 



Appendix B: Kinetic Energy 



We use a choice of the phase of the d-orbitals so that the starting kinetic energy in the model is a real operator. 
With this choice and the notation specified in Fig.([T|), the kinetic energy operator is 



K.E. = X] ~9~ (4iPii^ + d'liPi'tx - d^Piix - dl^Pi4x + 4iP'iy + 42P'2y - dliPiiy ~ AzP^^y 



(Bl) 



+ U 



(p\ixPay - PayP^i^ + pIaxPhv ~ p\iyP^i^) + 



h.c 



r 




FIG. 4: The five areas in which a unit-cell is divided by 
connecting Cu-O and O-O links. Sum of directed Current 
operators on the links to form closed loops define five fiux- 
operators in the areas marked by Roman letters. These can 
be further combined to form 5 fiux patterns with the point 
group symmetry of the square lattice as in fig.©. 



Appendix C: Gauge invariant combination of the 
vector potentials 

The complex hopping inatrix elements on the links 
in the unit-cell are equivalent to vector potentials living 
on the links of the unit-cell. Flux operators are formed 
by sum of the phase difference (or vector-potentials) in 
closed loops of links which form independent areas in 
each unit-cell. 

There are 12 Cu-O and 0-0 links shown in Fig.([T|) but 
8 of these are shared by the adjoining unit-cells, so that 
there are only 8 links per unit-cell. But the links form 
only five areas per unit-cells as shown in Fig. ((4]). This 
is consistent with the fact that there arc three sites per 



unit-cell in which independent gauge transformation can 
be made to obtain 5 gauge invariant flux-operators per 
unit-cell. 

The five closed loops in any unit-cell can be combined 
to form five new coinbinations which with the point- 
group symmetries of the lattice. The five combinations 
are pictorially represented in figlS] The five combinations 




FIG. 5: The five gauge invariant combinations of link vari- 
ables that respect the symmetry of the lattice. Li^i has the 
symmetry of an overall fiux in an unit cell. Li, 2 has zero to- 
tal fiux in an unit cell. The corresponding operator has the 
symmetry of an angular momentum. Uniform ordering of Li, 3 
results in the S/ phase. Li,4 and Li, 5 are the operators that 
condense to give the Qn phase. 



of /lus-operators L^ „ are, 



-L^J.t' 



i.x'^ —y'^ 



^i^x' 



i-'i,y' 



I \d\^piix - d\^pi4,,^, + d\2Piix - dl^pi4x - dliPiiy + d^^pi2y - d\^piiy + d^ap^yj 
I {d\^p^ix - d\^piix + dl^Piix ~ dj^piix - 4iPiiy + d^Pi^y - d^Ptiy + dlaPt^y) 

« {p\xP2y - plyPlx+pLpiy - piyPlxj + h.C 

I f -dJiPila; - dl^Pi4r, + d^aPilx + d\r^p^i^ + d\-^Pily + d\2Pi2y - dj^Pily - d\r^Pi2y 

* \p\xP-2y + plyPix + p\xPly + p\yPlxj + h.C 

i {dl2Pilx + dl^Pi4x + dl^Piiy + dl2Pi2y] + C.C 

* [PlxP2y - pLpiv) + h.C 

i [dl-i^Piix + dl^pi4r, - dl^pay - dl^pi2y] + h.C 

« [PlyPix - plyPlxj + h.C 



h.C 
h.C 



(CI) 



As will be evident from Fig. ([5]), Lg has the identity 
representation of the flux-operator, while Lg- in common 
parlance may be called the "extended s-wave" represen- 
tation. Lix' and Ly' are the operators that have the sym- 
metry of the 6// phasei. The phases produced by their 
condensation have magneto- electric symmetry, describ- 
able by time-reversal odd polar-vectors pointing in the 
x',y' = (x±y)/"\/2 directions respectively. 'Lx'2_yi2 has 
the symmetry of the 8/ phase described earlier—. 



Appendix D: Derivation of the Interaction 
Hamiltonian in terms of flux operators 



In this appendix, we give details of the re-cxprcssion of 
the Cu-0 and 0-0 interaction Hamiltonian iJ„„ given by 
Eq.(IT]) in terms of current operators using the operator 
identity of Eq. ^. The purpose of doing this is to 
derive the relevant part of the interaction Hamiltonian 
Hnn in terms of the irreducible combinations of the flux- 
operators of Eq. ([3]) . The unimportant one-electron terms 
in Eq.Q arc ignored and only the spin-diagonal parts 
are kept. The spin-diagonal part of the interaction H„n 
across the 12 links in Fig.([T]) 



here: 



^"« = -^^[|o..,.p + ia 



(Dl) 



K 



-^ Z^ \^hi,xy 



Here i sums the four Cu-sites per unit-cell and the op- 
erators O, given after Eq. ^ are written down again 



O, 



i,£,x 



> ''■'^te c,Pi,t,x,a 






h.C, etc, (D2) 



O. 



i,i,x,y 



E 



wt,i,y,aPi.^.=:.<y + '^■c-, etc. (D3) 



To simplify notation, we use a slightly modified label- 
ing scheme in this appendix. In a given unit cell there 
are four triangles with one of their vertices being a Cu 
site and one square with its vertices being the four Oxy- 
gen atoms. The subscript x and y refer to the Px and 
Py orbital, that coixibined with the Copper site, labeled 
by £, form the triangle. Explicitly the triangles are 

{l,Plx,Ply}, {2,pix,P2y}, {3,P4x,P2y} and {i, P4x, Ply} . 

The flux in the triangles labelled L = /, ...IV in Fig.(|3]) 
with a Cu-sitc at £ = 1, ..4 is, 

fi,I = OiS,x ~ Oi^i^y + Oi^i^xy (D4) 

fiJI = Oi^2,x + Oi^2,y + Oi^2,xy 

fiJII = —Oi^3,x + Oi,3^y + Oi^z^xy 

fi.IV = —Oi^4,x — Oi,4^y + Oi^4^xy 

A clockwise choice of currents around the loop is chosen 
to define / to be positive. The convention is shown in 

figEl 

Similarly, the fifth flux, fiy operator (See Fig. ([5])) is 



fi,v — 2^ fi,i 



20, 



,xy 



(D5) 



Here the sum of the currents on the Cu-0 links cancels 
the current in the 0-0 links such that there is no fiux in 
the corner triangles. 
For the triangle L = 1 

'Vpd 
4 ^ 

|2\ 



^(IO.mP 



I^M.yP) 



\{0^. 






O, 



(D6) 







FIG. 6: Convention used for defining tfie triangles and tfie 
link operators. 



We now note that the first term in Eq. (|D6p has finite sum 
of currents at the Cu-sites and can form closed loops only 
by adding similar terms from neighboring cells. These 
cannot give rise to either q = order or long wavelength 
fluctuations and are not considered further. Now we add 
{Vpd/4:)\Oi,e,x.y\'^ to the second term in Ea. (|D6p and sub- 
tract it from the similar term with coefficient Vrm in 



Eg. ljDip so that we can write 



pp 



H, 



L=l 



\{Oi,i^x — Oi.l,y ~ Ois,x,y)\ ) ■ 

{2Vpp - Vpd) 2 
-A l<-^»,l,a:y| 



Note using Eq. (jD4p that the first term in Eq. (jD7|) is equal 
to ffi^^i- This exercise can be repeated for L = 2,3 
and 4. For each triangle we get three terms: 1) ff^^, 
2) term analogous to the second term in Eq lD4l and 3) 
{{2Vpp — Vpd) / '^)\Oi_i,.xy\^ ■ We can sum over t in the sec- 
ond term to produce one combination which (see Eq. (jD5[) 
is |/i,yp. The other three can be removed by gauge 
transformations at the 3 sites in each unit-cell. 
Next note that 



Li,£c' — fi.i — fi.III 
^i,y' = fi,II — fi,IV 



(D8) 



L,,./2_„/2 — fi.I + fi.III — fiJI — fi, IV 



Li,s - 77 / , /i 



L=IJV 
i,L 



L=I,V 



We also note that 



Ho 



i,t,x.y ^i,s 



(D9) 



as well. 

Using the above identities, we can write the gauge in- 
variant part of Hnn as 



Hnn = -{^)[\U,x'\' + \Uy'\ 

+ (^)|Li,x'2-y'2| -t- |Lj_s| 

Vpi 
8 



This is in the desired form. 
Let us also define an operator 



(DIO) 



Li, 2 — 2^ fi, 



L=I,..,IV 



A term 



Hke — — 



16 ' 



(Dll) 



(D12) 



is also present in the interactions. We are concerned in 
this paper with the fluctuations of the observed phase 
which is realized by the local condensation of the collec- 
tive parts of Lix' and Tuyi into four possible domains in any 
unit-cell. In the quantum-fluctuation regime, the impor- 
tant fluctuations are between these four conflgurations in 
any unit-cell. As we show below, such fluctuations are 
caused by the operator L^ so that Hke acts as the ki- 
netic energy operator rotating the conflgurations of L^;/ 
and Ly' . The relevant terms in the above must therefore 
include Hke beside those involving the L^^', hy' operators. 
This completes the derivation of Eq.Q of Sec. II. 

The operator L^^s corresponds as noted to a net flux in 
the unit-cell. For it to order, its expectation value must 
reverse between neighboring cells giving rise to flux pat- 
terns with broken translational symmetry. Such patterns 
are not observed in experiments. What about the uni- 
form ordering of L^ s? This would correspond to a net 
flux in the sample or macroscopic boundary currents. By 
a general theorem, such ordering is impossible because 
long-wavelength variations of L^^^ have all the symme- 
tries of magnetic field produced by a vector potential, 
which cannot acquire mass. However from the same con- 
siderations, fluctuations of such operators at long wave- 
lengths have the properties of photons and therefore may 
be quite important in the pseudogap phase. The prop- 
agator of such fluctuations and its coupling to fermions 
and the consequences of this coupling will be presented 
in the near future. 



Appendix E: Properties of Li,z 

In this appendix we show that the operator L^ defined 
on the lattice is a generator for rotations. To do so let 
us consider the four states represented pictorially in figH] 
and write down the one particle wavcfunction which has 



10 



the same expectation value of the current operator as 
that in the collective states Lj in a unit cell. Define a 
set of basis operators centered at each of the four copper 
sites in a cell i labelled by their center, \ipii)^^ , |V'i2)^^, 

V'^^" = ^{V2di3 + e"l'p,4^ - e"''ip,2y) 

i'ii' = i^i^dii + e"''P^ix - e"''ip.ay) 

(El) 

In terms of these operators, the wavefunction of the 
electron corresponding to, say, the state 16* = 225°) = 
\-x - y) of Lj is 

\~x ~y)^\ {i\~ + iT + i^r + ^r)^ |0) (E2) 

Similar representation of the other four-states in terms 
of fermion operators may be written down. The state 
1^=135°) = \~x + y) is obtained by rotating clock- 
wise by 90". The corresponding rotation operator is 



U (-7r/2) = e-*L. W2) = I _ jL, - l' 



Using 



[d^d^ + d'^di,ii2] = J31, 
and [ji2, J23] — iji3 we can explicitly verify that 



U (-7r/2) , i (V^f- + iq- + V3"- + ^4—)^ 



(E3) 



(E4) 



(E5) 



(Vi + + '02 ^ + V^3 ^ + ^4 



-+\t 



where 



4^-+ = i(\/2d,i + e^Va. + e^Vay) 



li^d,. 



i2 + e"^pii^ 



o«0 



i^. 



-+ 
3 



e ^PiAx 



o'0 



iPi2y) 
iPi2y) 



fpii^ = ■^{V2di4 + e'-'^'piix + e'"^ipiiy). 



(E6) 



The exercise can be repeated to show that the operator 
indeed rotates the states of L clockwise by 90°. The 
lattice fermion operator L^ is the generator of rotation 
and corresponds to the angular momentum operator in 
the continuum limit. 



Appendix F: Correlations of the angular 
momentum operator 

The coupling of the fermions to the fluctuations of 
the order parameter are proportional to the angular 
momentum operator L^. The pairing interaction that 
leads to superconducting instability involve the (LziJjzj) 
correlation function. To analyze the connection with 
the {exp{i9{r, t)) exp(i0(r', t'))) correlation derived in Aji 
and Varma^ we first must identify the appropriate repre- 
sentation of the field exp {t9{r' ,t')). In the path integral 
formulation the correlation function computed is 



a. 



(^ 



ixg- 



^') 



(Fl) 



where 0^^ is the classical phase at site i at time fi. In this 
calculation the standard procedure of slicing time into in- 
finitesimal segments is followed. The commutation prop- 
erties are accounted for by appropriately defining matrix 
elements for infinitesimal time evolution as dictated by 
the Hamiltonian. We now show that the correlation func- 
tion of the angular momentum operators is proportional 
to the same classical field correlations. 

To do so we note that the angular momentum in its 
eigenbasis \jn) is given by 



hz = / m I'm) {ml 



(F2) 



We can now look at the operation on a state 16*) at site i 
and time r, 



IfzirlOir) = ^mirlm^r) {m\0ir) 

rrii-T- 

= ^77i,;,|m,,)e*«- (F3) 



The correlation of the angular momentum operator in 
the theta basis is 



E 



zir'-'zJT' 



(F4) 



-ZTTli-rO, 



^+™'^-'V' {niirlmjr') 



mi-T- ,vri- 



Since the (exp(i(?(r, i)) exp(i0(r',i'))) correlations are lo- 
cal the same is assumed to be true for the angular mo- 
mentum. This is justified because in the critical regime 
the spin waves are not the critical modes and the overlap 
of the angular momentum eigenstates {mir\ itt-jt') which 
are spatially separated falls off exponentially. At the 
same site (mi,-! rriiT') = Smi^m ,■ Having obtained the 
representation of the angular momentum correlation in 
the 9 basis, we can compute the correlation function. 
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FIG. 7: The spectral function of the autocrrelation as a func- 
tion of , shown. For m's larger than 1 the spectral weight 
shifts to higher frequency leading to short range correlation 
in time. For ni=l the spectrum leads to a power law decay. 



faster in time so that the dominant contribution to L^L^ 
correlations is precisely of the form in cqnfT51 We refer 
to the derivation in ref.— for further details but quote 
the main result. The correlations are given by 



(exp(l771r T^(r, r)) CXp(l77lr' t'S{v' , t'))) CX Sy. r' CXp V '^ J 

(F6) 
and 



F{r 



-2TTTm- 



1 -cos(a;„ (r- r')) 



log {uj„ 



(F7) 
where w„ is the matsubara frequency and Tc is the short 
time cutoff in the theory. For m ~ I, the long time cor- 
relations decay as |r — r'| . For m > 1, the correlations 
decay faster as the weight is shifted to higher and higher 
frequencies. One can infer this from the spectral function 
associated with this autocorrelator which is given by^i 



Cf 



L'^i^i-'z 



E ^" 



(F5) 



Ijr' t 



T+irn- 



The 



= 1 and 



IjT-l 



1 contribution to the angular 
momentum correlations in the sum in egn lFSI is equal 
to (exp(z6'(r,i))exp(z0(r',i'))). Thus the leading term is 
local in space and power law in time. Higher order in 
m correlations are also local in space but decay much 



A, 



coT'- 



sinh(#) 



r(^ 



27rT 



2(m=-2)/2^r (ml) r ( 



rn^ + l ' 
2 , 



(F8) 



where cq is a constant. To show the shift in weight we 
plot the spectral function for m = 1, 2, 3 and 4 in fig jT] 
As m becomes larger the frequency dependence at low 
frequencies is superlinear. The increasing weight at larger 
frequencies implies that the these correlations will decay 
faster that the m = 1 contribution. 
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